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ABSTRACT 

We present hydrodynamic simulations of molecular outflows driven by jets with 
a long period of precession, motivated by observations of arc-like features and S- 
symmetry in outflows associated with young stars. We simulate images of not only H2 
vibrational and CO rotational emission lines, but also of atomic emission. The density 
cross section displays a jaw-like cavity, independent of precession rate. In molecular 
hydrogen, however, we find ordered chains of bow shocks and meandering streamers 
which contrast with the chaotic structure produced by jets in rapid precession. A 
feature particularly dominant in atomic emission is a stagnant point in the flow that 
remains near the inlet and alters shape and brightness as the jet skims by. Under 
the present conditions, slow jet precession yields a relatively high fraction of mass 
accelerated to high speeds, as also attested to in simulated CO line profiles. Many 
outfiow structures, characterised by HH222 (continuous ribbon), HH240 (asymmetric 
chains of bow shocks) and RNO 43N (protruding cavities) , are probably related to the 
slow precession model. 

Key words: hydrodynamics - shock waves - ISM: clouds - ISM: jets and outflows 
- ISM: molecules 



1 INTRODUCTION 

Collimat ed outflows are int e Rrally connected with st ar for- 
mation jRicher et al.l I2OO0I: iReipurth fc Bailvl I2001II . The 

outflowing mass is distributed over a large volume and can 
far exceed the suspected mass of the growing star. This sug- 
gests that momentum is generated deep within the potential 
well of the young star and that, after extraction, is trans- 
mitted into the ambient gas. Jets of molecules are often as- 
sociated with the youngest protostars. However, one view 
suggests that these jets are highly ballistic and that the low 
density of the surroundings does not efficiently soak up the 
jet momentum. As the protostars age and the accretion rate 
declines, their molecular jets become less dense and more 
atomic which may increase the efficiency. Alternatively, the 
jet thrust might be sprayed into a large volume through pre- 
cession of the direction of release. The consequences of such 
a scenario are explored here. 

Observations of protostellar jets frequently provide ev- 
idence for precession. Some examples with precession an- 
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gles i n the ran ge 5° - 20° are Cep A fGueth (?^&ulloteau 
1999[) : Cep E llEisloffel etallllQQd); HH 31 (IGomez et al. 
19971): HH 34 |Ballv fc Devind Il994^ : and HH 333 
T^cefc Goodm an^O^) . Larger precession angles a re sug- 
geste d for L 1228 (BaUv ct al. 199h). PV Cep (Gomez et alJ 
Il997^ and IRAS 20126-H4104 lshepherd et aljboOOl) . 



We present here a sequence of three dimensional hydro- 
dynamic simulations. They were executed with a modified 
version of ZEUS-3D which solves a detailed cooling function 
and a reduced chemical network (see . We drive a molec- 
ular jet into a uniform molecular environment. Previously, 
we have examined the physical and observational structure 
which results for jets with (i) a range of den sities relative to 
the ambient density jRosen fc Smithll2004lJ) . (ii) significant 
s usta ined changes to the inlet mass flux (|r osen fc SmithI 
[2003^ and (iii) a range of precession angles bu t all with 
a fast precession period jRosen fc SmithI l2004a^ ■ Fast jet 
precession, however, does not produce the S-symmetry as 
observed in the above examples. Instead, the jet impact de- 
scribes an expanding ring which fragments into a multitude 
of shocks. In addition, a stagnation region appears once the 
impact no longer encounters the precession axis, as exam- 
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ined by iLimI (I200J) for an atomic model. The ambient gas 
trapped in this region may lead to observable features. 

Here, we investigate the consequences of jets under- 
going slow precession, the precession speed being defined 
in terms of the appropriate dynamical time (see ij^Jl. We 
wish to identify the global structures associated with a 
range of commonly observed emission lines. Curved out- 
flows of various types are indeed observed. The very young 
protostar VLA 1623 is associated with a slender tube-like 
CO outflow which meanders through the p Ophiuchi cloud 
iDent et al.l lll995l) . A few bipolar outflows show C-type or 
mirror symmetry rat her than S-type or po int symmetry (e.g. 
the S233IR outflow. iBeuther et al.ll2002^ . Such flows could 
result from twin jets associated with one protostar within a 
close binary: the space velocities of both jets at each instant 
are the vector additions of fixed components in the pr oto- 
star's frame and the variable orbital velocity (Frank o t alJ 
Il998l: iMasciadri fc Raeal l2002ll . There are even some the- 
ories that tie the precession of jets to the warping of the 
disk iLaill20o3) . The star formation process, which must be 
able to generate a large number of binaries, thus includes 
potential mechanisms for varying the direction of ejection. 

Besides identifying outflow properties, there are two is- 
sues we try to resolve. Firstly, strong H2 1-0 S(l) emission 
is often measured in the near-infrared at 2.12 /im. Taking 
this to be ten percent of the total H2 emission, which is 
taken to be ten per cent of the total radiation, then the jet 
power must be at least 100 times stronger. This can just be 
reconciled with current theories for young. Class proto- 
stars which may drive jets of high power in comparison to 
their bolometric luminosity. However, taking near-infrared 
extinction into account as well as the fact that ballistic jets 
tend to be inefficient radiators, the jet powers become un- 
comfortably high. 

Another motivating measure concerns the distribution 
of outflowing mass with radial velocity. All of our previous 
three-dimensional simulations predict shallow distributions: 
a high fraction of gas is accelerated to high speed. In con- 
trast, there are numerous outflows for which steep relation- 
ships are deduced from CO line proflles. We summarise these 
measurements of the slope in the mass- velocity and intensity 
distributions in 5 15. II We test here if these problems with the 
efficiency and velocity distribution might both be solved by 
slowing the rate of precession. 

The sim ulations are executed with a modified version 
of ZEUS-3D JSmith fc Rosenll2003ll which tracks the molec- 
ular fraction, calculates a quite detailed cooling function, 
and includes molecular dissociation and reformation. The 
cooling and chemistry are so lved implicitly. Details and 
tests have been presen ted by ISmith fc RoserJ ll2003ll and 
iRosen fc Sraitiil (l2004d) . 



Table 1. Major parameters altered between the simulations dis- 
cussed (upper set) and the initial/boundary values for the ambi- 
ent medium and jet (lower list) which were not varied. 



Run 


PIO 


QIO 


RIO 


Precession period 


50 


400 


400 


Precession angle 


10° 


10° 


10° 


Pulsation period 


60 


60 


60 


Pulse amplitude 


30% 


30% 


0.01% 


Jet speed (axis) 




100 km s-i 




Jet speed (perimeter) 




70 km s^i 




Jet density 




2.32 10-19 gm cm-3 




Jet temperature 




100 K 




Jet radius 




1.7 10^5 cm 




Ambient density 




2.32 10-2" gm cm-3 




Ambient temperature 




10 K 




n(H2)/nH 




0.5 




Init. specific heat ratio 




1.4286 
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2 ZEUS-3D WITH MOLECULAR COOLING 
AND CHEMISTRY 

2.1 Hydrocode properties 

We have modified the ZEUS-3D code, which usually up- 
dates variables with an explicit scheme, to include a semi- 
implicit method to calculate the time-dependent molecular 
and atomic hydrogen fractions. We have added a limited 



Figure 1. Midplane cross sections of density. Each cross section 
is scaled to the same maximum and minimum and the scale is 
displayed below the cross sections. The complete computational 
domain in the two axes is displayed. The vertical axis is the z-axis 
in each panel. 
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Figure 2. Axial cross sections of density for simulations PIO (top row), QIO (middle row) and RIO (bottom row). Each row of snapshots 
is taken from the time that is displayed in the corresponding panel in Fig.Q The axial position in Rj is given in the upper left of each 
panel. Each cross section is scaled to the same maximum and minimum and the scale is displayed below the cross sections. The complete 
computational domain in the two axes is displayed. The distance between adjacent tick marks is 1 X 10^^ cm. The vertical axis is the 
y-axis and the horizontal axis is the z-axis, and so the view is from the jet inlet boundary looking down the x-axis (i.e., toward the 
+x-direction). 



equilibrium C and O chemistry to calculate the CO, OH 
and H2O abundances. Equilibrium CO and H2O is a rea- 
sonable estimate at high density, consistent with the low 
spatial resolution of the cooling layers behind shock waves 
jRosen fc SmithI ^004b) . The details of the many compo- 
nents of the cooling function and the li mited chemistry 
netwo rk are discussed in the appendices of ISmith fc RosenI 
i2003h . In summary, we include cooling through rotational 
and vibrational transitions of H2, CO, and H2O, H2 dis- 
sociative cooling and reformation heating, gas-grain cool- 
ing/heating, and a time-independent atomic cooling func- 
tion that includes non-equilibrium effects. We assume that 
the gas includes an additional 10% by number of helium 
atoms, so the mean atomic mass is 2.32 x 10~^* g. The 
dust temperature is fixed at 20 K. 

The cooling and chemistry are poorly resolved as a re- 
sult of the demands placed by the need for three dimen- 
sional simulations. The problem is located within the ra- 
diative cooling zones behind shock waves. Since fast shocks 
with speed exceeding ~ 25 kms~^ dissociate molecules, both 
atomic and molecular cooling zones arise in the present sim- 
ulations. In fast shocks, the atomic cooling length is con- 
siderably shorter than the molecular cooling length and so 
provides the most stringent scale to resolve. However, this 
length is of the order of just 10^^ - 10^* cm, given the initial 
range in density. Therefore, we make no attempt to resolve 
the atomic cooling zone and note that peak temperature 



obtained in the simulations will be considerably less than 
the physical values associated with the shock speed. For- 
tunately, as utilised in p revious studies (|Stone fc Normanl 
liggdlsuttnei" et al .119971) ■ shock conservation laws are accu- 
rately maintained despite the lack of resolution. This is used 
to argue that the dynamics are correctly predicted although 
the atomic emission line fluxes are incorrect (the images of 
total atomic cooling presented here should be considered as 
providing only the general locations of shocked atomic gas) . 



The molecular cooling zones are significantly longer: ~ 
10^" /n cm. While shocks into the lower-density ambient gas 
are resolved, those in the jet are generally not. Hence, we 
restrict our predictions to the properties of the global out- 
flow rather than individual shocks. This is particularly rel- 
evant to the chemistry. By choosing equilibrium H2O and 
CO chemistry, the influence of the cooling zone on the entire 
flow is minimised. However, a non-equilibrium treatment of 
H2 is necessary and, unless care is taken, this could lead to 
misleading results. Tests carried out by us at various levels 
have shown that the resolution achieved here is adequate to 
describe the global properties but that much more detailed 
fine structure will appear at higher resolution iSuttner et alJ 
I1997I: iPavlovski et alJl2002^■ Other potent ial causes of inac- 
curacy were discussed bv lRosen fc SmithI J200l l2004d) . 
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2.2 Chosen simulations 

In a previous work, we described a simulation of a molecular 
jet with a precession angle of 10° and a fast precession rate, 
which we designated PIO. In this paper, we discuss PIO and 
two additional simulations with a 10° precession angle, but 
with a much slower precession rate (see Table 0. One sim- 
ulation is identical to PIO save for the precession rate, and 
is designated here as QIO. The other simulation is identical 
to QIO with the exception of the amplitude of the jet pul- 
sation, which is reduced from ±30% to ±0.01% (essentially 
a constant initial jet velocity). We designate this simulation 
RIO. 

We modify the jet velocity in PIO and QIO with a 30% 
pulsation and a 60 year period (see Table0. In addition, all 
simulations have a radial shear that reduces the jet velocity 
at the jet edge to 0.7 Vj. The precession period is 50 yr in 
PIO, but 400 yr in QIO and RIO. The former is relatively 
slow compared to the dynamical time Rj/vj = 5.4 yr, but 
fast relative to the flow evolution time of order 500 yr (see 
below). The slower rate is of the same order as the flow 
evolution and also our a;- axis crossing time. 

We define the precession rate as fast if it rotates faster 
than the outflow expands. Then, the jet contributes to ad- 
vancing an entire annulus into the ambient medium, gen- 
erating a ring structure or annulus of fixed width (twice 
the jet radius). To calculate the expansion time, we first 
note that the density of the precessing jet decreases to 
^ Pj/{1 + {D sin 0)/rj) at a distance D along the x axis. This 
is effectively a one-dimensional expansion since the width of 
the annulus remains constant as its distance from the axis 
about which the jet is precessing increases. Momentum bal- 
ance with the external medium of uniform density pa then 
yields 



1 + {D sin e) jr. 



(1) 



on equating the ram pressures, where U is the average 
advance speed of the fast precessing jet into the ambient 
medium. If we now consider the expansion to be braked 
when U is reduced, say, to 0.5 Vj, then the expansion time, 
tE ~ D/U, is given by 



Pi - Pa 



2r, 



pa 



(2) 



If the precession is slow then we expect a meandering 
stream to have time to develop, rather than a ring. For the 
parameters chosen here, ts = 559 yr. Therefore, our cho- 
sen 'slow' precession rate still displays the effects due to 
precession whereas a 'gradual' precession, with a period far 
exceeding would approximate to a non-precessing struc- 
ture. 

The computational volume is simulated with a uni- 
form staggered grid divided into zones each of which spans 
2 X 10^* cm in each direction. For all three of the 10° 
precession simulations, we use a grid of 480 x 205 x 205 
zones. Since the jets have large precession angles, the even- 
tual time-averaged momentum axis for the fiow (here, the 
a;-axis) is never the jet flow direction. The initial jet radius, 
Rj, is 1.7 X 10^^ cm and is therefore resolved by 8.5 zones. 
The boundary zones satisfy outflow conditions with the ex- 
ception of the zones on the inner a;-boundary that allow jet 
inflow. 



We assume conditions that we suspect are relevant to 
Class protostellar jets. In all cases, the jet is initialised 
with a hydrogenic nucleon density of 10^ cm~^. We set the 
jet-to-ambient density ratio to 10, the initial jet temperature 
to 100 K and the ratio of jet-to-ambient thermal pressure to 
100. The nominal, mean axial jet speed is Vj = 100 km s~^. 
This yields an initial Mach number with respect to (the uni- 
form) ambient medium of 140. 



3 PHYSICAL PROPERTIES 
3.1 Internal structure 

Midplane slices of density are displayed in Fig.Qat a sim- 
ulation stage where side projections of the outflows almost 
spanned the entire computational grid. In fact, at the times 
shown in Fig.Q both QIO and RIO have jet material as far 
from the inlet as a: ~ 8.5 x 10^^ cm (50 Rj). Most striking is 
the existence of a jaw-shaped structure in the slow precess- 
ing cases QIO and RIO as well as the fast precessing case. 
Also, slow precession generally yields lower density cavities 
(darker in Fig.0 since swept up material is not quickly re- 
placed by the jet. The midplane density slices reveal the 
shocks from the pulsing in QIO, and relative smoothness of 
the non-pulsed flow in RIO. 

We display the cross-sectional structure normal to the 
jet axis for the same three simulations in Fig.|21 For PIO, this 
shows a complete annulus mid-way down the jet which frag- 
ments near the leading edge. For QIO and PIO, the slow 
precessing results in a helical flow pattern from +z (the 
right side of each panel) to —z through +y (the top of each 
panel) and back again through —y. This is consistent with 
the time- variation of the precession at the origin. A compari- 
son between the QIO and RIO jets reveals that the instability 
seen in the "islands" of jet material at large x are primarily 
caused by the pulsing. Speciflcally, the islands seen in QIO 
at x/Rj = 31 and 39 are unifled into an arc-like bow struc- 
ture in RIO. However, as presented below, the integrated H2 
emission portends future fragmentation in the RIO flow. 

Quite high transverse velocities are encountered in the 
slow precessing cases. A plot of velocity vectors and density 
contours (Fig.|2J for RIO shows that outward facing velocity 
vectors are associated with the helical bow shock far from 
the inlet. Closer to the inlet (in the first couple of panels), 
the precessing jet is evident. In addition, there is a large 
transverse velocity (~ 50 km s~^) in the third panel, where 
the bow shock has nearly spun around on itself. 



3.2 Speed and size 

We have previously analysed the fast precessing jet simula- 
tions to determine their capability of prod ucing giant out- 
flows of size ~ 1 pc iRosen fc Smithll2004aD . We employed a 
model to fit the mean speed of advance and to determine the 
effective drag coefficient corresponding to the aerodynamic 
properties and atmospheric resistance. This term, Cd, can 
also be parameterised as a change in area across the bow 
shock in the usual ram pressure approximation for the Mach 
disk advance speed, as used in §|5| We deduced low values of 
Cd in the range 0.18-0.27 for a range of density ratios and 
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Figure 3. Transverse velocity vectors and mass density contours 
in simulation RIO at t = 350 yr. The value labelled x is the ax- 
ial position and Vmax is the maximum transverse velocity in the 
panel (in km s~^). The velocity vectors are displayed for every 
sixth zone and the maximum vector length in any panel is nor- 
malised to the maximum velocity in that panel. Density contours 
are at log p = -24.0, -23.0, -22.0, -21.0, -20.0, -19.0, and -18.0. 
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Figure 4. Advance speed and estimated time to 1 pc for precess- 
ing molecular jets. The panels on the left show the time-averaged 
(over 20 years) advance speed of the bow shock, which are indi- 
cated by the triangles connected with the dashed line. The ad- 
vance speeds have been fitted with a model that includes the 
momentum drag from the interaction of the jet with the ambient 
medium (Rosen & Smith 2004a). The inferred drag coefficient, 
Cjj, is listed in the upper right of these panels. The models are 
indicated by the solid lines, the most accurate includes the slowly- 
damped sinusoidal variation caused by the precession while the 
monotonically decreasing fit assumes that the jet momentum is 
evenly spread in the advancing disc/annulus. The panels on the 
right show the time required to reach 1 pc, according to the ex- 
trapolation of the advance speed. 



mass fluxes whereas we expect a value of unity if the ram 
pressure of the jet is modified by its expansion only. 

For the fast precessing cases, the drag coefficient 
increases with precessing angle, from 0.20 for 5° half- 
precession angle to 0.45 for the 20° case. Within this model, 
we found that the fast precessing jets could reach a dis- 
tance of 1 pc between 35,000 and 85,000 yr, showing ev- 
idence for deceleration despite the small drag coefficients. 
We concluded that outflows reach large sizes in times com- 
parable to the Class protostellar stage. 

Here, the evolutions of the linear extents are shown in 
Fig-Bl Also shown are the best fit models and integrated 
times to 1 pc. We find that the drag coefficient varies ac- 
cording to the physical structure. The advancing bow shocks 
of run QIO result in a low drag coefficient of 0.22. In con- 
trast, the continuous ribbon produced by run RIO leads to 
the highest drag (0.45). This is consistent with expectations 
from fluid dynamics since the 3D aerodynamic bow shapes 
permit fluid to flow easily past the obstacles, whereas the 
ribbon presents a more obtrusive object. The time required 
to reach 1 pc according to the extrapolated model remains 
within the 40,000-60,000 yr range. 



4 SIMULATED MOLECULAR LINE EMISSION 
4.1 Integrated line emission 

We estimate the associated radiated energy rate for a few 
molecular and atomic emission lines using the zonal values 



for mass density, temperature and molecular fraction. We 
integrate the emission for selected molecular lines at spe- 
cific times with the assumption that the outflow and fore- 
ground cloud are optically thin; examples of these integra- 
tions are shown in Figs.lS HTTI The H2 line emission is based 
on a non-L TE approximation to the vibrational populations 
jPraine et a l. 1983) and the CO r otatio nal emission is cal- 
culated according to Mc Kee et al.l(ll982^ . The only plot pre- 
sented for run PIO (Fig.|^, which was discussed previously, 
includes new integrations in "atomic" emission as well as in 
molecular hydrogen emission for comparison. 

The atomic emission is de rived from the cooling func - 
tion used in the code (based on lSutherland fc Dopitalll993l) . 
unmodified to select a specific emission line. As with all 
of these integrations, they provide a guide to the location 
and extent of emission and relative brightness of regions of 
the fiow rather than accurate predictions. In run PIO, the 
brightest atomic emission is spatially correlated with rela- 
tively bright regions of molecular hydrogen emission, a result 
that is not surprising since both are related to the cooling of 
post-shock material. However, the brightest atomic emission 
occurs in localised small knots and filaments whereas the H2 
emission is relatively more diffuse. 

The images of atomic emission prominently display in- 
ner cones associated with a stagnation point (the knot at 
3 X 10^^ cm and some of the emission forward of that). 
Also, there is no extended atomic emission from the lead- 
ing edge of the jet, where there is significant H2 emission. 
This suggests that the atomic emission is produced in strong 
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Figure 5. Two orthogonal views of simulation PIO at t = 450 yr 
in one molecular emission line and one general "atomic" emission 
line. In the top two panels, the integration is along the j/-axis, 
and in the bottom two, integration is along the z-ajds. In order 
to produce a more complete picture of atomic emission, the im- 
age has been rescaled so that the maximum luminosity in the 
appropriate panels is 3 orders of magnitude lower than the actual 
maximum luminosity (which is very localised). The luminosities 
are in erg s~^ from bins with roughly the same size as the 3D 
zones used in the simulations (i.e. 2 X 10^* cm). 
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Figure 6. Integrated emission along the y-axis in simulation QIO 
at t = 350 yr for four molecular emission lines and one general 
atomic emission line. Further details are as in Fig.|Sl 



shocks directly associated with the jet impact rather than 
the predominantly weak oblique shocks associated with the 
advancing aerodynamic bow-shaped features. 

The simulated images of molecular hydrogen and 
atomic emission reveal the shock structure. In simulation 
QIO (Figs. |S| and [7|, the pulsing jet at the inlet leads to 
many internal shocks, while in RIO (Figs.|H]and|31l the con- 
stant speed yields a shock in the form of a single sinuous 
ribbon. 

This ribbon contains an apparent "x" at a:: ~ 2.5 x 10^® 
cm in the y-integrated images (Fig.|SJl. It is most prominent 
in the atomic emission from simulation RIO but is also bright 
in all emission lines shown. While more difficult to perceive, 
a similar feature is also distinguishable in the plots for run 
QIO. In the images generated by z-integration (Fig.jnj, this 
"x" is a more complicated structure, suggestive of a disk of 
emitting material seen edge-on. Such a disk is confirmed by 
viewing the jet from various angles about the mean momen- 
tum (x) axis, which we display in Fig. llOl Moreover, this disk 
can be seen in in the second panel of the RIO density cross 
sections in Fig.|5|and contours of the third panel in Fig.|3| 

From these images, we hypothesise that the "x" is re- 



lated to a stagnation disk that at the time of this image 
is being squeezed from both sides, since the precession rate 
is just sufficient for more than a half-turn to exist at this 
a;-position. In addition, the third panel in the velocity vec- 
tor plot has a very large maximum transverse velocity (~ 
50 km s ~^ compared with 17.5 km s~^ transverse velocity 
from the nominal jet speed and the input angle). It seems 
plausible that this large velocity is caused by the release of 
accumulated material from this stagnation disk. To investi- 
gate further, we have also executed isothermal simulations 
of a fast dense jet which show evidence for large tangential 
velocities near the jet-crossing point. 

The "x" in atomic emission appears similar to that in 
the HH222 (seen in a plate of the IIH34 superjet complex, 
iBallv fc Devinelll994ri . although the the process responsible 
is likely different in the real source, which is many times the 
size of what has been simulated here. 

We present the temporal variation of molecular hydro- 
gen emission from run QIO in Fig. llll As opposed to the fast 
precessing cases, where at early times we saw an elliptical 
ring form, here we see a strong asymmetry about the mean 
momentum axis. However, at early times this axis may not 
be obvious. Additionally, the bow shocks evolve in a compli- 
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Figure 7. Integrated emission along the z-axis from simulation 
QIO at t = 350 yr in four molecular emission lines and one general 
atomic emission line. Further details are as in Fig.lSl 
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Figure 8. Integrated emission along y-axis from simulation RIO 
at t = 400 yr in four molecular emission lines and one atomic 
emission line. Further details are as in Fig.l^ 



cated and fascinating fashion, with the strongest emission 
moving forward through the bow as the entire emitting re- 
gion moves forward. For times just after t = 300 yr, there 
is a rapid change in the direction of the jet, which could be 
(mis-)interpreted as a dramatic event at the jet origin. Of 
course, in this case we know that only a smooth variation is 
responsible. The sequence of b ow shocks at late times bears a 
striking similarity with L1634 JO'ConneU et al.l2004^ . which 
is an order of magnitude larger than our simulations. 

The vibrational excitation is measured by the H2 2- 
1/1-0 line ratio. The excitation distributions reveal the most 
intense shocks, and its brightest regions generally coincide 
with bright regions of the 1-0 line itself. Noteworthy is that 
the region associated with the "x" feature in RIO has quite 
a high H2 line ratio. 

The moving CO emission displays an outline of the 
dense jet plus swept-up ambient material. The CO emis- 
sion lines, particularly the 0-0 R(l), emphasise the sinu- 
soidal projection of the helical flow pattern. The ribbon is 
traced back to the source. Note that the individual advanc- 
ing bow shocks in run QIO sweep out separate protruding 
CO cavities (Fig.|SJ. Finally, the CO R(5) images (i.e. in the 
rotational emission line with a frequency of 691 GHz, from 



116 K above the ground state) again appear somewhere be- 
tween the H2 emission and the CO R(l), as expected. Both 
the cavity structure and internal bullets are evident. 

4.2 Total line emission and efficiency 

The evolution of the total molecular line emission from 
molecular outflows is of relevance to the global modelling 
of star forma tion ( e.g., the Unifica tion Scheme, |Smith 20(3 
IStankdIioOol: and IYu et ai]|2000f) . In contrast with previ- 
ous pulsating and some of the fast precessing simulations, 
the total molecular intensities for molecular hydrogen and 
CO monotonically increase over time for the new simula- 
tions presente d here. They cover the same range as run PIO 
(see Fig. 10 in lRosen fc Smit'lill2004al) . The non-pulsing case 
RIO increases in the most smooth fashion of the three sim- 
ulations discussed here. The periods of flat or decreasing 
luminosity that exist in QIO coincide with the appearance 
of a new knot near the inlet. Eventually this knot brightens 
sufliciently that the total flux increases. 

We have also calculated the ratio of the total emission 
compared to a time-averaged kinetic energy flux (mechanical 
power) from the jet source for a sample of our molecular jet 
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at t = 400 yr in four molecular emission lines and one atomic 
emission line. Further details are as in Fig.lSl 
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Figure 10. Close-up rotation of views about the axis in simu- 
lation RIO. By rotating about the jet axis of the "x" feature in 
Fig.|S] we see a more varied appearance. 



simulations (see Fig. 11 from lRosen fc Smitbll2004'ah . Typi- 
cally, the H2 1-0 emission line had an efficiency of around 
0.1-1.0% of the mean mechanical energy, with an increase 
in this range over the lifetime of the simulation. For CO the 
range varied from 10"'^ at early times to roughly 10~^ at late 
times {t ~ 500 yr), with a more monotonia variation with 
time between the two extrema than in the H2 case. We have 
also computed the evolution of these ratios for the two new 
simulations and find values close to that of Simulation PIO. 
In general, both the variation over time in the mechanical 
power and in the H2 to mechanical power ratio (by 0.2 dex, 
usually) is smaller in the slow precessing cases than in PIO. 
Also, the ratio is typically smaller (0.1 dex) in the pulsed 
case (QIO) than non-pulsed case (RIO). For the CO ratio, 
QIO is smaller or equal to the ratio in PIO, and the RIO 
CO ratio is typically larger or equal to that of PIO, but over 
time both asymptotically approach the evolution of PIO. 

We have also calculated the the instantaneous radiative 
loss from each part of our cooling function at 50 yr inter- 
vals for the new simulations. These losses are then compared 
with the time-averaged mechanical luminosity (power) input 
at the inlet up to the point in time that the losses are com- 
puted. We find that at t = 350 yr, the total losses are roughly 
30% of the mean power, with two-thirds of that being from 



a combination of H2 ro- vibrational and atomic cooling. For 
the pulsed case (QIO), the atomic emission alone is about 
5% of the mean power at this time, but for the non-pulsed 
case (RIO), it's fraction is only 0.5%. This is explained by 
the many internal bow shocks, which are sites of heating and 
dissociation for the initially molecular gas. The evolution of 
the H2 1-0 S(l) emission efficiency from 0.1% to 1% in these 
simulations is associated with a nearly constant 4% ratio of 
emission from this line to the total H2 ro-vibrational cooling 
rate. 

Other major contributors to the cooling function in- 
clude (values in parentheses are the fractions of each cool- 
ing component of the mean power for the jet near the end 
of the simulation): H2 dissociation (2-3%), H2O rotational 
transitions due to collisions with H and H2 (~ 5%), and CO 
rotational transitions due to collisions with H and H2 (2%). 
The individual component ratios of luminosity to mean jet 
power nearly always increase monotonically with time, as 
does the total, which increases by about 5% of the mean 
power in every 50 yr interval. 
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Figure 11. Evolution of tlie 1-0 S(l) molecular liydrogen emis- 
sion for simulation QIO. Integrations of H2 1—0 emission along 
the 2-axis (the y-axis is the vertical dimension in all the panels). 



The simulations in iRosen fc Sniithl J2004bl) possess flux- 
velocity distrib utions of H2 that do not closel y resemble the 
observations of'Salas fc Cruz-Gonzale3 i2002l) . with a strong 
increase in luminosity at small velocities. This is not dis- 
turbing since the 1-0 S(l) emission properties are very sen- 
sitive to the physics within the radiative shocks. In hydro- 
dynamic shocks, the molecules radiate efficiently only after 
being strongly accelerated and compressed. This suggests 
that the shock viscosity is instead provided by ambipolar 
diffusion in the molecular sections. 

A xisymmetric simulations simila r to those in 
iDowne s & Rav ( 1991) were analysed by iDownes fc Cabril 
i2003r) . They obtained a low-velocity CO 7 of 1.5, which 
is in the range of commonly observed values. This may 
have been the result of smaller jet densities and longer 
time coverage than we have used. Indeed, ISmith et alJ 
(1997) also found that the 7 valu e incre ases with time in 
3D simulations. IDownes fc Caj3ri^_l|200 ^ also attempt t o 
reproduce the H9 results of ISala^gCru^-Gonzaled i2002f) . 
with some success but use an approximation (LTE) that is 
more appropriate for densities larger than in most of their 
flow. 



5 VELOCITY DISTRIBUTION OF MASS AND 
MOLECULAR LINE EMISSION 

5.1 Background 

There have been many recent measurements of the slope 
in the mass-velocity distributions, interpreted from CO line 
profile data of protostellar outflows. Once the slow turbu- 
lent motions of nearby molecu lar clouds are meticulously 
subtracted fe.g.. lYu et aljlioool) . there is frequently a quite 



shallow distribution (with 7 in dm/dv 



occasionally 



as low as 1) at low radial velocities. A steeper slope ( ^ 6) 
is found at higher velocities, where the turnover occurs at 
Vbreak '~ 10 km s~^. Notc, however, that the majority of 
powerful ou tflows may posse s s very steep profiles even at 
low speeds JPavis et al.lll998l: iRidee fc Moorell200J) . Many 
of these outflows possess low coUimation, suggesting that 
the wide outflow sweeps up a high mass to low speeds. 

Additionally, a similar analysis has been per- 
formed for H2 lin e profiles of some protostellar outflows 
JSalas fc Cruz-Gon zalcz 200j). The equivalent 7 is defined 
as 7m = dL{l — S{l))/dv oc v The authors measure 
a range of 7^ between 1.8 and 2.6 above a break velocity 
which lies between 2-17 km s~^ (the flux is flat or slightly 
rising below this radial velocity). We wish to test here if jets 
undergoing wide-angle long-period precession could deceler- 
ate sufficiently to match the observations of large j,n at high 
speeds. 

Previous simulations of ballistic jets with small pre- 
cession angles yielded quite shallow values for the CO 7 
of 1.2-1.6 (Sm ith ct al. 1992) or even shallower (0.9-1.3, 
iRosen fc SmithI l2004bl) . Non-ballistic jets, with equal or 
lower density than the am bient medium, also gen erated out- 
flows with low 7 values l|Rosen fc Smith l2004bl) . The fast 
precessing jets in iRosen fc SmithI J2004al) yielded a some- 
what larger range of slopes for the small velocity region, 
(0.4-1.6). The higher 7 values are associated with combina- 
tions of the following conditions: jets viewed near the plane 
of the sky, jets with a large precession angle, or jets which 
are initially atomic. 



5.2 Mass- and intensity- velocity relations 

We have performed a similar analysis for the velocity distri- 
butions of mass, CO R(l) luminosity and H2 1-0 luminosity 
for each of the simulations presented here. Power law slopes 
7, 7co and yH2 are determined. We list the results for a sam- 
ple of viewing angles, defined as the angle of the a;-axis out 
of the plane of the sky toward the observer, in Table |5| We 
also display the data for one viewing angle (15°) in Fig. 1121 
Note that both red and blue shifted line wings are displayed 
in this figure. A summary of our results follows. 

• The precession period has a small influence on the mass 
distribution as measured by the power-law index 7. For 
viewing angles > 15°, there is a slightly steeper slope in 
the mass relation in the slow-period precessing cases, which 
have slopes that are fairly close to each other (the non-pulsed 
case marginally increases the slope). 

• For the CO distribution, there is a stronger and oppo- 
site dependence of yco on the precession rate, with most 
viewing angles of the slow-precession cases yielding a slope 
0.5 smaller than the rather shallow 1.2 for PIO. 

• In simulations of fast precessing jets, we found that the 
slopes for the CO distribution exceed those of the mass 
distribution, occasionally quite significantly. This is sim- 
ilar to a result from t he low density jet simulations of 
IDownes fc Cabril i2003l) . This trend extends to the slowly 
precessing molecular jet simulations presented here, al- 
though the differences between the two slopes are small in 
QIO and RIO than in PIO. 

• At least for large viewing angles (> 45°), we confirm 
previous results that there is an inverse relationship be- 



tween 7 and the viewing angle 1 Smi^^^L|^9^jLe^^i] 
I2OOI', 'Downes fc Cabriti 120031 iRosen fc SmithI l2004hl and 



£,osen fc Smith 2004a). 

• Additionally, we have examined the evolution of the ve- 
locity distributions for each of the molecular jet simulations 
at two different programm e times. In contrast to the find- 
ings of lSmith et al.l I^1997^ ■ we find that the slope of these 
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Table 2. Mass Spectra Power-Law Dependences 



Sim. 


t(yr) 


type 


e 


range of log \v\ 


7 




type 


e 


range of log \v\ 


7 




PIO 


450 


mass 


15 


0.5-1.0/1.0-1.4 


1.02/1.87 


7/15 


CO 


15 


0.5-1.0/1.2-1.4 


1.20/1.96 


7/9 








30 


0.5-1.0/1.1-1.7 


0.52/1.68 


7/37 




30 


0.5-1.0 


1.22 


7 








45 


0.5-1.0/1.4-1.8 


0.63/1.65 


7/38 




45 


0.5-1.0 


1.22 


7 








60 


0.5-1.0/1.4-1.9 


0.57/1.45 


7/54 




60 


0.5-1.0 


1.21 


7 








90 


0.5-1.0/1.6-1.9 


0.46/1.37 


7/39 




90 


0.5-1.0 


1.21 


7 


QIO 


350 


mass 


15 


0.5-1.0 


0.26 


7 


CO 


15 


0.5-1.0 


0.49 


7 








30 


0.5-1.0 


0.56 


7 




30 


0.5-1.0 


0.76 


7 








45 


0.5-1.0 


0.63 


7 




45 


0.5-1.0 


0.80 


7 








60 


0.5-1.0 


0.59 


7 




60 


0.5-1.0 


0.74 


7 








90 


0.5-1.0 


0.53 


7 




90 


0.5-1.0 


0.65 


7 


RIO 


400 


mass 


15 


0.2-0.7/0.8-1.3 


0.35/1.12 


3/14 


CO 


15 


0.2-0.7/0.8-1.3 


0.83/1.35 


3/14 








30 


0.5-1.0 


0.63 


7 




30 


0.5-1.0 


0.71 


7 








45 


0.5-1.0 


0.72 


7 




45 


0.5-1.0 


0.79 


7 








60 


0.5-1.0 


0.65 


7 




60 


0.5-1.0 


0.76 


7 








90 


0.5-1.0 


0.56 


7 




90 


0.5-1.0 


0.70 


7 



^ N is the number of points in the velocity distribution, which was computed in 1 km s ^ bins, used to estimate the slope. 



distributions shallow slightly over the small length of time 
simulated here. 

• The velocity distributions of H2 emission for the slow 
processing cases contrast with the obscrvationally deter- 
mined ones, as with our previous simulations. Instead of 
the observed flat distribution at low velocities, we find a 
dramatic rise. 



and 



5.3 Discussion of position-velocity images 
channel maps 

In previous papers, we have presented position-velocity (P- 
V) diagrams in molecular emission lines from various sinm- 
lations. From simulations of slow precessing jets, we find no 
additional features and so omit any new figures here. 

However, there are some differences. The main one is 
that in the middle of the projected jet (x') axis, the H2 
diagram appears much less "turbulent", even in the QIO 
pulsed case, than in PIO and thus the fiow is much better 
focussed for the slower precession rate. In the RIO case, the 
middle portion is contained completely within 0-10 km s~^, 
although this is not true for the CO emission, which extends 
out to 20 km s~^. 

In RIO, rather than the sequence of Hubble Law like 
P-V regions (that are a consequence of the pulsing) we find 
an arc of bright material that starts at high velocity (35 km 
s~^) near the origin and swings down to 10 km s~^ near 
the middle of the plot. About half-way along (from about 
one-quarter the jet length in x' out to the maximum extent 
of the jet), the entire range of velocities out to the arc are 
filled in. In addition, the maximum velocity increases from 
the middle portion outward, to 25 km s~^. 

Analysis of previous velocity channel maps, primarily 
in CO, revealed a morphological progression as the viewed 
velocity is changed. High velocity maps show the isolated 
knotty structures associated with internal shocks along the 
jet and low velocity maps display the full extent of the 
shocked ambient medium. In the fast precessing cases, the 
knots in the high velocity maps are displaced from the mean 



-2 h 

° -3 t 

^ -5 h 
a> i 

o -6 r 

-7 h 



Q 

< -10 
o 
o 

f -12 

-Uh 

_ 3 
-4 

o 

^ r 

-7 |- 



-9 h 



0.0 



58 Run QIO, t = 350 yr 
+ Run RIO, t = 400 yr 



* ******** 
* 



+ ++* 

+ 

I 1 1 1 1 1 

^ SIS * si< it slisl?**±, 



+++ 




0.5 1.0 1.5 

log IvrodI (km/s) 



2.0 



Figure 12. Distribution of velocities into bins of mass and two 
molecular line luminosities for a viewing angle of 15°. The distri- 
butions of mass (top), CO luminosity (middle) and H2 1-0 emis- 
sion (bottom panel) arc shown. Each velocity bin is 1 km s~^ 
wide. The designation for the data presented in each panel is 
shown in the top panel. Each run may be represented twice within 
each panel, with data from both positive and negative radial ve- 
locities included. Naturally, the lower valued data is for the pos- 
itive radial velocities (which could contribute to a "red" lobe, 
while the other is from a "blue" lobe). 
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momentum axis, frequently in an alternating pattern. In the 
slow precessing, pulsing jet case (QIO), we see a sequence of 
alternating knots, but the spatial frequency of the oscillation 
has decreased. In both QIO and RIO, the velocity channel 
maps indicate that the bright arc within the CO emission 
near the inlet is at high velocities, while much of the inte- 
grated CO emission (Figs.|Sland|Hl is at low velocities. 



6 SUMMARY OF RESULTS 

We have found that the dominant physical structure for all 
precessing jets is an inward facing cone. In other words, the 
jet clears out a similar volume of ambient gas. The jaw- 
shaped structure contrasts with the tube-like structure ex- 
pected if the jet were to be bent or deflected. 

A fast precessing jet generates an impact ring which 
fairly rapidly disrupts into many mini-bow shocks. Here, 
slowly precessing jets lead to helical flows that remain no- 
ticeably more stable for the duration of the simulation. 

A slow precessing and pulsating jet cuts through the 
ambient medium with the highest speed. We measure this 
here through an effective drag coefficient. The non-pulsating 
jet results in a drag coefficient about double that of the 
pulsating case. Therefore, bow shocks represent the most 
efficient means of advancing an outflow yet the least efficient 
at transferring momentum. 

The stagnation region at the apex of the inward fac- 
ing cone displays quite complex time-dependent behaviour 
in the slow precessing cases. The feature does move slowly 
forward and accumulates some material, which after some 
time can be sloughed off to the side open to the ambient 
medium. We believe that this is the cause for the high trans- 
verse velocities measured in these simulations (particularly 
the non-pulsed RIO case). The feature, which at times can 
resemble a disk in 3D, is seen in atomic and high energy 
molecular emission as an "x" feature from certain viewing 
angles. 

The images of molecular emission show quite dramatic 
evolution. In the pulsed case, the H2 emission shows a train 
of asymmetric wedge-shaped bow shocks. The sequence of 
bow shocks can appear over time to quickly change direc- 
tion suggesting a violent event at the driving source, yet 
we are simulating a smoothly varying jet flow. The distri- 
bution and structure of the bow shocks appear very sim- 
ilar t o the observed HH24 0A-D system in the L 1634 out- 
flow dO'ConneU et all2004^ . Note, however, that the HH 240 
bows have advanced a distance of ~ 0.4 pc from the driv- 
ing source. Many other groups of bo ws hav e been uncovered 
(e.g. HH195 and HH197 5 (iEisloffellll2000D ). 

On the other hand, shock emission in the form 
of ribbons is rare ly encountered. Examples arc HH 222 
jBallv et alJ l2002tl. HH 103-105 (Eisloffel 2000) and 
IRAS 23151+5912 Iweigeh et"aiT2005h . In fact, HH 222 also 
possesses an "x" feature, resembling that predicted in our 
simulations. 

The CO cavity patterns found here are largely feature- 
less. In the pulsed case, however, distinct shells may outline 
the wakes of separated bow shocks which advance directly 
into the ambient cloud. To observe this property requires 
resolution and sensitivity which may soon be generally at- 
tainable. At present, only one well-known precessing out- 



flow, RN043, is k nown to possesses this striking property 
jBence et al.lll99dl . 

Many bipolar outflows contain a high fraction of cool 
gas moving at radial velocities near to that of the embedded 
cloud. This is derived from CO intensity-velocity relation- 
ships which are approximated by a high negative power-law 
index. In contrast, low power-law indices are associated with 
these simulations. Hence, the solution to the discrepancy be- 
tween our results and observations lies elsewhere. It might 
be reconciled by running the simulations for longer peri- 
ods of time, which should be easier to perform as faster 
and larger computing facilities become available. Despite 
the lack of resolution some interesting outstanding issues, 
we have shown that the precession rate can dramatically in- 
fluence the appearance of precessing molecular protostellar 
jets. 

The simulations ignore the effects of the magnetic field. 
Simulations under ideal magnetohydrodynamic conditions 
may show significant differences. For example, the mag- 
netic field may stabilise the compressed layer produced by 
the advancing shock wave. In this case, the fragmentation 
into multiple bow shocks could be inhibited. In addition, 
there is mounting evidence that ambipolar diffusion is im- 
portant i n protostellar ou tflows iDavi j I2OO2I : ISmith et alJ 
2003: Gia nnini eEaOEooi). Ion-neutral friction allows the 
magnetic fleld to efficiently cushion the shock, spreading and 
combining the compression and cooling zones. This would 
have two major effects. Firstly, molecules are not so effec- 
tively destroyed until shock speeds approaching 50 kms~^ 
are reached. Hence, the atomic emission will be even more 
tightly distributed than here. Secondly, molecular emission 
is generated from a wider range of shocks. This will tend to 
increase the length of the wings of bow shocks and the flanks 
of cavities. Hence, MHD simulations which are directly com- 
parable to those presented in this series should provide more 
insight. 
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